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Abstract. We derive the expression for the thermal conductivity k, in the low-temperature limit T — > in 
d-wave superconductors, taking into account the presence of competing orders such as spin-density wave, 
is-pairing, etc. . The expression is used for analyzing recent experimental data in La2- a; Sr I Cu04. Our 
analysis strongly suggests that competing orders can be responsible for anomalies in behavior of thermal 
conductivity observed in those experiments. 

PACS. 74.25.Fy Transport properties (electric and thermal conductivity, thermoelectric effects, etc.) - 
74.72.Dn La-based cuprates - 74.72.-h Cuprate superconductors (high-Tc and insulating parent com- 
pounds) 



The existence of four nodal points in d-wave super- 
conductors provides rich and, sometimes, controllable dy- 
namics of quasiparticle excitations at zero temperature. 
In particular, the expressions for electrical, thermal, and 
spin conductivity simplify considerably in the universal- 
limit us — > 0, T — > PH]. It is noticeable that the role of 
the thermal conductivity k is special: while vertex and/or 
Fermi-liquid corrections modify the bare, "universal" , val- 
ues of both electric and spin conductivities, the univer- 
sal value of the thermal conductivity is not influenced by 
them 2]. It is: 



explicitly, a nonzero Jo is crucial both for Eqs.Q and 

But what will happen if those quasiparticles become 
gapped? One may think that in that case both N(0) and 
Krj are zero. However, as will be shown in this paper, they 
both are finite even in that case, if the impurity scattering 
rate is non-zero. In fact, it will be shown that they are: 
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where vp is a Fermi velocity, Va is a gap velocity, and 
ks is the Boltzmann constant (we use units with fi 
c 1 . . The basis for such a remarkably simple expression 
is that there is a finite density of states N(Q) of gapless 
quasiparticles down to zero energy [2|3]: 
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where To = r(uj — > 0), with r(to) an impurity scattering 
rate, and po = tJttvfva/ a is an ultraviolet momentum 
cutoff (a is a lattice constant) 0. Note that expression 
CJ itself is valid in the so-called "dirty" limit, T -C T . 
Therefore, although this expression does not contain iy 
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where m a quasiparticle gap. The noticeable point is that, 
for all values of the gap up to m ~ To, the suppression 
of both thermal conductivity and quasiparticle density is 

mild: K /K m) and N(0)/N m (0) are of order one. However, 
the suppression in thermal conductivity rapidly becomes 
strong as m crosses this threshold. The second noticeable 
point is that, as we will discuss below, the gap m plays 
here a universal role and may represent different compet- 
ing orders in d-wave superconductors, such as as spin den- 
sity wave, charge density wave, is-pairing, etc. . Although 
their dynamics are different, expressions © and for 

k^ 1 " 1 and N rn are the same. This happens because, first, 
all those gaps m correspond to different types of " masses" 
in the Dirac equation describing nodal quasiparticle exci- 
tations, and, secondly, unlike electric and spin conductiv- 
ities, the thermal conductivity fc m ' ) is blind with respect 
for Theoretical ^° quantum numbers distinguishing those masses. 

The expression corresponds to the dirty limit 

when T <C Tq. In d-wave superconductors, Tq can be as 
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large as of order 1 K or even 10 K, and fc^™ 1 ' can be an im- 
portant measurable characteristic there. Recently, two ex- 
perimental groups have observed an anomalous behavior 
in the thermal conductivity in underdoped L^-zSr^CuC^ 
(Refs. and Refs. |MZj)- One of the most interesting 
observations of experiment 6 is that at very low tempera- 
tures the value of the thermal conductivity in underdoped 
La2- 2: Sr a: Cu04 (LSCO)is less than the absolute minimum 
kmin/T — 2fcf 3 /3 of expression (JTJ for ko/T, correspond- 
ing to the isotropic case with Vp — va- This puzzle can 
be naturally explained by utilizing the modified expres- 
sion © with a nonzero m describing a competing order 
in the superconducting phase. We will discuss this and 
other results of experiments 0|EllSllZj below. 

At subkelvin temperatures relevant to the low-T heat 
conduction experiments, we will use the continuum, low- 
energy, description for the nodal quasiparticles in the d- 
wave state. At each node, the quasiparticles are described 
by a two-component Nambu field. It will be convenient, 
following Ref. (Hj, to utilize four-component fields, by com- 
bining Nambu fields corresponding to the nodes within 
each of the two diagonal pairs. Thus we have two four- 
component Dirac fields. The corresponding representation 
for three Dirac matrices is 



7o = 0i ® I, 7i 



-ia 2 ® cr 3 , 72 = ia 2 ® <J\ , 



(5) 



where a% are the Pauli matrices and while the first factor 
in the tensor product acts in the subspace of the nodes in 
a diagonal pair, the second factor acts on indices inside a 
Nambu field. The matrices satisfy the algebra {7^,71/} = 
20/xz,, 6V = (1, -1, -1), (J*,v = 0, 1, 2. 

We will consider quasiparticle gaps with the matrix 
structure O; = (I, 275, 73, 7375). Here the matrices 73 and 
75, anticommuting with matrices j u , are 



73 = «0- 2 ® 02, 75 =03 



I. 



(6) 



Go(iuj n — E(iu) n ), k). As usual, retarded Green's function 
is obtained by analytically continuing Green's function G, 
G r (oj, k) — G{iuj n — ► u) + ie, k), and the impurity scatter- 
ing rate is defined as r(ui) = -ImS R (w). At low temper- 
atures we take i~o = r(ui — > 0). The size of Jo depends 
on the impurity density rii mp as well as on the scattering 
phase shift S. Solving the Schwinger-Dyson equation for 
the self-energy in the self-consistent i-matrix approxima- 
tion, one can find that in the unitary limit (5 = 7r/2) the 
equation determining Jo for a nonzero m, has the form 



r 2 = n 2 v F VAr 
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where Nf is the number of four-component Dirac fields 
and r — riimp/irpQ with po the normal state density of 
states. Since va ~ Ao, the magnitude of the supercon- 
ducting gap, the scattering rate Zo is proportional to 
V A r ~ A n lmp . 

The longitudinal dc thermal conductivity is calculated 
by means of the Kubo formula. In the bubble approxima- 
tion, following the standard procedure, it can be expressed 
through the quasiparticle spectral function A(u), k) as fol- 
lows 
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Here the spectral function is given by the discontinuity of 
the fermion Green's function 

A(u, k) = -— \G R (u + ie, k) - G A (u - ie. fell . (10) 

Z7TI 



Then, for each of the two four-component Dirac fields, the 
bare Matsubara Green's function can be written as 



G {iui n ,k) 



1 
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Therefore, different gaps rrij correspond to different types 
of Dirac masses. As was pointed out in Refs. |9lllllll()l 
IHj, these gaps represent different competing orders in low 
energy limit. In particular, the mass mi, with 0\ = I, 
describes the (incommensurate) cos spin-density-wave 
(SDW), and the mass m%, with O2 — 275, describes sin 
SDW. The masses 7713 and 7714, with O3 — 73 and O4 = 
7375 j correspond to the id^-pairing and the is-pairing, 
respectively. One can also consider a gap corresponding to 
the charge-density- wave (CDW). In that case, one should 
introduce a Dirac mass term mixing the four-component 
Dirac fields corresponding to the two different diagonal 
pairs of the nodes. For simplicity, we will not consider it 
in this letter. 

The scattering on impurities can be taken into ac- 
count by introducing a Matsubara self-energy E(ioj n ), so 
that the dressed Green's function becomes G{iui n ,k) = 



With Green's function at hand, we can calculate A(u>, k). 
For example, for the gap proportional to the unit Dirac 
matrix, it has the form (m = mi) 
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where E(k) = \/v 2 F k\ + v 2 A k\ + m 2 is the quasiparticle 
energy. Substituting the last expression in Eq. @ and tak- 
ing the limit T — ► 0, we arrive at 
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i.e., we derived expression Q for the thermal conductivity 
(in which Nf — 2). The result for three other gaps, m 2 , 
TO3, and TO4, introduced above, is the same. 
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With the spectral function the density of states 
(per spin) 

N m (u) = I [ j^tr [ 7o A(cj, fc)] (13) 



2 J (2tt) 2 



is easily calculated 



N m (u) 



N 



f 



To In 



rn In • 



y/r* + (u;-my 



Po 



_ 1 u 2 — m 



7T 

— + tan 



r 2 

1 o 



2Ur 



(14) 



It yields expression Q for the density of states with zero 
energy. Therefore, in the presence of impurities, the quasi- 
particle band survives even for a finite m. 1 The physical 
reason for this is the formation of impurity bound states 
inside the gap Overlap between these states leads 
to impurity band supporting the quasiparticle heat (and 
electric) current. 

The observation of a residual linear in T term in the 
thermal conductivity in cuprates (YBa2Cu3 07_a (YBCO) 
[El, Bi 2 Sr 2 CaCu 2 8 (Bi-2212) ,16; as well as LSCO ©) is 
usually interpreted as a direct consequence of nodes in the 
gap. However, as it follows from Ea. H12|) . a subdominant 
order parameter, leading to a gap for nodal quasiparticles, 
does not exclude such a linear term in the thermal con- 
ductivity, although the latter does not have a universal 
form anymore. 2 

Thus we conclude that nonperturbative dynamics, re- 
sponsible for the creation of competing orders in the su- 
percritical phase, can violate the universality in the ther- 
mal conductivity in the low temperature limit T — > 0. 
Recent experiments indicate that the existence of such 
competing orders is quite possible [T7j . Several theoretical 
models have been proposed to describe this phenomenon 
(for a review, see Ref. QB]). As we will now discuss, using 
the expression for the thermal conductivity derived above, 
this phenomenon can be relevant for understanding recent 
experiments in L^^-zSr^CuC^ 0]|S]|()]|7| . 

The measurements of the thermal conductivity in 
LSCO at low temperature 4,5 6 7 showed the following 
characteristic features: 

a) At subkelvin temperatures, the value of k/T de- 
creases with x 0JE1- At temperature as low as 40 mK, the 



1 Note that in the absence of impurities [To = 0], we would 
get N m (u>)\r =o = (Nf /2-kvfva)\oj\0(uj 2 — m 2 ), i.e., in that 
case, the mass m would lead to a gap in the density of states. 

2 Although the fact that opening of a gap for nodal quasipar- 
ticles leads to changes in Kg™' ^ s natural from physical view- 
point, there has been a controversy concerning this point in the 
literature. For example, in the recent paper the authors 
claim that in the limit T — > the universal expression for the 
thermal conductivity, Eq. survives even for gapped quasi- 
particles. Expression 112H derived above clearly shows that this 
is not the case. 



value of k/T in some underdoped samples is either less 
than the absolute minimum K m i n /T — 2k 2 B /2> of expres- 
sion (JJJ (for x — 0.06) or quite close to it (for x — 0.07 and 
x = 0.09) UJ- On the other hand, this anomalous behav- 
ior in the thermal conductivity disappears in overdoped 
samples (x = 0.17 and x = 0.20) 0. 

b) The evolution of k/T across optimum doping is 
smooth UK]. 

c) The thermal conductivity is sensitive to magnetic 
field. While in overdoped samples it increases with mag- 
netic field, in underdoped samples the thermal conduc- 
tivity decreases with increasing magnetic field [HHZ]- The 
authors of Refs. 017] describe this as a field-induced ther- 
mal metal-to-insulator transition. 

d) Although remaining smooth, the evolution of k/T 
across optimum doping becomes visibly faster with in- 
creasing magnetic field [S]- 

The results of item a) can be easily understood if one 
assumes that there exists a competing order, described by 
the Dirac mass m, in the superconducting phase of un- 
derdoped LSCO. Then an appropriate value of m in ex- 
pression @ will provide the necessary suppression of the 
thermal conductivity. The fact that such an anomalous 
behavior in k/T disappears with increasing x, in over- 
doped samples, can be understood if one assumes that 
the dynamical gap ("mass") m decreases with increasing 
x. As to this assumption, it is well known in quantum 
field theory that, indeed, an increase of the fermion den- 
sity often suppresses a dynamical Dirac mass. The reasons 
for that are simple. With increasing the fermion density, 
the screening effects become stronger and the quasipar- 
ticle interactions become weaker. In addition, at a suffi- 
ciently large quasiparticle density, the energy gain from 
creating a gap m in the quasiparticle spectrum will be 
surpassed by the energy loss of pushing up the energy of 
all states in the band above the gap. In the case of the 
model with Dirac fermions describing highly oriented py- 
rolytic graphite (HOPG) 191120) . this fact was explicitly 
shown in Ref. 20 . Although the present system is quite 
different from HOPG, that example supports plausibility 
of this assumption. 

It is tempting to speculate that the dynamical gap m 
disappears close to optimum doping (xq = 0.16 in LSCO). 
A smooth evolution of k/T across optimum doping then 
suggests that it could be a continuous phase transition 
with the scaling law of the form m ~ (x c — x) v in the scal- 
ing region with < [x c — x)/x c <C 1, where the critical 
value x c ~ xq. The critical index v = 1/2 would cor- 
respond to the mean-field phase transition. In that case, 
there would be a kink in expression © at the critical point 
x = x c . Indeed, since the thermal conductivity © de- 
pends on m 2 , and there is a linear in m 2 term as m 2 — > 0, 
its derivative with respect to x will have a finite disconti- 
nuity at x = x c for v = 1/2. In the case of a non-mean-ficld 
continuum phase transition, with v > 1/2, the evolution 
of k/T across x c ~ xq would be smoother. 

This picture, with appropriate modifications, can sur- 
vive in the presence of a magnetic field. In particular, the 
fact that in overdoped samples k increases with magnetic 



4 



V.P. Gusynin, V. A. Miransky: Thermal conductivity and competing orders in d-wave superconductors 



field as y/~H |5][7], implies that the dynamics in a magnetic 
field in overdoped samples is apparently conventional. In- 
deed, the y/H behavior is well described by semiclassical 
models This seems to suggest that there is no gap m 
(competing order) in overdoped samples. 

The situation is different in underdoped samples. The 
magnetic field enhances the suppression in k observed in 
the same samples at zero field (item c) above). Moreover, 
the evolution k/T across optimum doping becomes visi- 
bly faster with increasing H (item d)). This suggests that 
magnetic field plays here the role of a catalyst, enhancing 
the gap m. For sufficiently large values of m, the suppres- 
sion in k will be so large that a sample effectively becomes 
a thermal insulator as was observed in experiments |UCj- 

Microscopic dynamics responsible for creating compet- 
ing orders can be quite sophisticated ^H]. This issue is out- 
side the scope of this letter. Here we will only comment on 
the role of a magnetic field as a catalyst in generating the 
gap m. In non-superconducting systems, it is well known 
that a magnetic field is indeed a strong catalyst in gener- 
ating gaps (masses) for Dirac fermions [22] • In particular, 
this effect was studied in the model describing HOPG |191 
I2()j . It is clear, however, that the dynamics in the vortex 
phase of d-wave superconductors is very different and the 
question about the relevance of a magnetic field for gen- 
erating (or enhancing) a quasiparticles gap there is still 
open. For example, while the authors of papers |23U24U25| 
believe that such a role for a magnetic field in that phase 
is plausible, the analysis of the authors of Ref. 2.6. indi- 
cates that the magnetic field can actually supress id xy and 
is gaps in a d-wave state. 

In this paper, we will use a heuristic approach and 
demonstrate that the experimental data in Refs. 017] can 
be qualitatively understood if one requires a gap that is 
generated below a critical doping and increases with a 
magnetic field. To make this point to be transparent, we 
are looking for an ansatz for the gap m( JJ, x) which would 
be as simple as possible. We assume that a) the phase 
transition at the critical doping x — 0.16 is the mean-field 
(or nearly mean-field) one, and b) the gap increases as \TS. 
with the magnetic field (such a scale covariant dependence 
of m on H was first considered in d-wave superconductors 
in Ref. [2H])- This leads us to the ansatz: 

m{H, x) = (l- x/0.l6) 1/2 61(0.16 - x)(m + bE H ), (15) 

where 9 is the step function, Eh — hvp/2R 
= (hvp/2)y/eH/hc is a characteristic energy scale in the 
presence of a magnetic field in the vortex state (2J? is the 
average distance between vortices) , and tuq and b are free 
parameters. Taking vp — 2.5- 10 7 cm/s for LSCO cuprates 
[TJETJ, we find E H = 38K-^/H(T) where the field H(T) is 
taken in Teslas. The constant b is of order 1 (for numerical 
calculations we take b = 2.2). As to the parameter mo 
that determines the gap for H — 0, it can be found from 
the ratio k/k q = 2/3 (i.e., k/T ~ 12^WK _2 cm _1 ) for 
x = 0.06, JJ = and T — > 0, as reported in Ref.®. Then, 
taking k/kq = Kq 1 /ko, with k™ from equation ©, we get 
mg = aJo where the constant a ~ 0.9. 



Let us now calculate the thermal conductivity by us- 
ing ansatz l|15|) for m(H,x). The impurity bandwidth for 
LSCO is estimated to be r ~ 25K - 30K pTl2Tj which is 
two orders of magnitude larger than for very clean 
YBa2Cu306.99 samples. While in clean YBa2Cu30e.99 the 
scattering of quasiparticles from vortices must be taken 
into account, one can neglect the dependence of the width 
Jo on the magnetic field (at least for not very high fields) 
in the case of rather dirty LSCO. On the other hand, the 
presence of a circulating supercurrent around vortices in 
the vortex state can be taken into account in the semiclas- 
sical approach by making the Doppler shift in quasiparti- 
cle energies, u) — ► to — v s (r)k, |2H] (v s (r) is the superfluid 
velocity at a position r which depends on the form of 
vortices distribution). In this case, the local thermal con- 
ductivity n(r) has to be averaged over the unit cell of the 
vortex lattice [21] , 



k(H,T) 



d 2 m{r) = / deV(e)K(e,T), 



where 



(16) 



(17) 



is the vortex distribution, and A — ttR 2 is the area of the 
vortex unit cell. We use the Gaussian distribution func- 
tion V(e) = (1/*/%Eh) exp[-e 2 /E H ] which is believed to 
be the most suitable distribution in the presence of high 
disorder jHO]. Thus we need to calculate 
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d 2 k 
(18) 



(compare with Eq.@). 

Taking the limit T — > in the last equation, we arrive 
at the following expression: 



k(JJ,0) _ 1 
k 2 



deP{e) 
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e 2 - to 2 + r 2 
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— tan 
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where we normalized the thermal conductivity on the uni- 
versal value Eq. . 

In Figs lll2l we present the ratio k(H,0)/k o calculated 
as a function of the magnetic field JJ (FigQ and the dop- 
ing x (Fig[2]). The form of these dependences is quite sim- 
ilar to that of experimental data presented in Fig. 2 of Ref. 
[7J and in Fig. 4 of Ref. 0, respectively. 

At small values of the doping, x — 0.06 and x = 
0.13 (low curves in FigQ, the thermal conductivity de- 
creases with increasing field as a result of increasing the 
gap m(JJ). For supercritical values of the doping (x = 
0.17; 0.2 - upper curves in FigQJ the field dependence is 
approximately \fH. This behavior is in accordance with 
the increase in quasiparticle population due to the Volovik 
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Fig. 1. k(H)/T [normalized to the value k(0)/T] versus 
H at T = and for the doping with the values x = 
0.06,0.13,0.17,0.2. The impurity width is To = 25K (for the 
upper curve (x = 0.2) 7~o = 30K). 



effect that is valid even for gapped quasiparticles |31| when 
the vortex scattering is neglected. 

FigEl shows the dependence of k on the doping for 
two different values of the magnetic field. One can see 
the suppression of k in the underdoped regime as a result 
of the presence of the magnetic-field-induced gap. Note 
that both curves grow fast near the critical doping x c = 
0.16 where the gap disappears. It is also noticeable that 
this growth is much faster for the H = 13 T curve than 
that for the H = 1 T curve. These facts agree with the 
experimental data :5 discussed in item d) above. 

Although the present analysis is based on the particu- 
lar ansatz (|15|l for m(H, x), one can expect that the main 
characteristics in the behavior of the thermal conductiv- 
ity will retain qualitatively the same for a wide class of 
gaps m{H, x) sharing the features that they are generated 
below a critical doping and increase with a magnetic field. 

In conclusion, we derived the expression for the ther- 
mal conductivity in d-wave superconductors in the pres- 
ence of competing orders. The derived expression @ for 

Kq" 1 ^ /T is simple and transparent. We also analyzed the 
dependence of the thermal conductivity on a magnetic 
field and a doping in the vortex state. Our results strongly 
suggest that the presence of competing orders can be cru- 
cial for understanding recent experiments in LSCO fJllHl 
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Fig. 2. Doping dependence of the T — thermal conductivity 
k [normalized to the universal value Ko] for two values of the 
magnetic field H — IT (dotted curve) and H = 13T (solid 
curve) and Jo = 25K . 
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